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ABSTRACT 

We compute exactly the frequency of Lense-Thirring precession for point masses in 
the Kerr metric, for arbitrary black hole mass and specific angular momentum. We 
show that this frequency, for point masses at or close to the innermost stable orbit, 
and for holes with moderate to extreme rotation, is less than, but comparable to the 
rotation frequency. Thus, if the quasi-periodic oscillations observed in the modulation 
of the X-ray flux from some black holes candidates, BHCs, are due to Lense-Thirring 
precession of orbiting material, we predict that a separate, distinct QPO ought to be 
observed in each object. 
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1 INTRODUCTION 

The large effective area, very high time-resolution and excel- 
lent telemetry of the Rossi X-ray Timing Explorer (RXTE) 
have made possible the discovery of Quasi Periodic Oscilla- 
tions (QPOs) in the range ~ 100 — 1200 Hz from a variety of 
accreting collapsed objects, weakly magnetic neutron stars 
(see van der klis 1998 for a review) and, more surprisingly, in 
black hole candidates (BHCs, Morgan et al. 1997; Remillard 
et al. 1997). It has been recently suggested ( |Cui et al. 199£ ) 
that these QPOs in BHCs arise through Lense-Thirring (LT, 
1918) precession of matter from the accretion disks. 

Since motion of a point mass in a Kerr metric allows an 
exact treatment, and a detailed comparison of Keplerian and 
Lense-Thirring frequencies has not been explicitly carried 
out in the literature up to now, it seems worthwhile to derive 
these quantities for an arbitrary black hole mass and specific 
angular momentum. This allows us to clarify the meaning of 
the precession frequency, about which some confusion seems 
to be present in the literature. This is the aim of this Letter. 
In the last section we shall also discuss the problems which 
are raised by this computation, regarding the interpretation 
of Cui et al. (1998). 



2 BOUND ORBITS IN THE KERR METRIC 

In what follows we will consider a test particle of unit mass 
in motion inside a Ke rr spacetime. The metr ic, in Boyer- 
Lindquist coordinates (Boyer & Lindquist 1967) and in units 



G = c = 1 is 
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M and a are, respectively, the mass and the specific angular 
momentum of the black hole. 



As Carter (1968) first demonstrated, the equation of 



motion can be separated, and the resulting equations be- 
come: 
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The dot denotes differentiation with respect to the proper 
time r; signs in (Q) and (^|) can be chosen independently. 
E, L and Q are the three constants of the particle motion: 
E and L are, respectively, the energy and the angular mo- 
mentum in the azimutal direction as seen by an observer at 
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rest at infinity; Q is relat ed to Carter's constan t of m otion 
(see e.g. Chandrasekhar (1983) and de Felice ( 198Cl| )) and 
characterizes the 9 mo tion. 

As Wilkins (1972) showed, bound motion is possible 



only if E < 1 and Q > 0; moreover, for given Q and L 
and \E\ < 1, there may be at most one region of binding. 
Analysis of the 9 effective potential shows that every orbit 
either remains in the equatorial plane (Q = 0), or crosses 
it repeteadly (Q > 0). For every bound motion, introducing 
the angle-action variables, we can define the three funda- 
mental proper frequencies 



1 / T~r,p 



where T<j> :P , t$ )P and r r ,p are the proper time periods for 
4>, 9 and r motions respectively. Unlike the Newtonian case 
of particle motion around a spherically symmetric central 
object, where all orbits close and the three fundamental fre- 
quencies are equal, in the Kerr field (a 7^ 0) there is no 
degeneracy, i.e. 

v 4>,p 7^ v e,v 7^ v r,v 

The same is also true for coordinate frequencies vg and 
v r . 

Let us first consider a circular geodesic in the equatorial 
plane (9 — tt/2). We have, for the coordinate an gular veloc- 
ities measured by an observer static at infinity (Bardeen et 
al. 197s|) 
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the angular velocity deviates from its Keplerian value 
at small radii. The upper sign refers to prograde orbits and 
the lower to retrograde ones. If we slightly perturb a cir- 
cular orbit introducing velocity components in the r and 9 
directions, we can compute the coordinate frequencies of the 
small amplitude oscillations within the plane (the epicyclic 
frequency Q r ) and in the perpendicular directi on (t he ver- 
t ical f requency Qg) (Okazaki, Kato & Fukue (1987), Kato 
(1990), de Felice & Usseglio-Tomasset ( |l996|) , Perez et al. 
(1997)): 
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In the case of the Schwarzschild metric (a = 0), there is a 
partial degeneracy, as the vertical frequency coincides with 
the azimutal one. The epicyclic frequency, instead, is always 
lower than the other two, reac hing a maximum for r = 8M 
and going to zero at r = 6M (Okazaki et al. 1987). These 



qualitative behaviour of the epyciclic frequency is preserved 
in the Kerr field (a / 0), and is a key feature fo r the e xistence 
of trapped diskoseismic g- modes (Perez et al (1997)). 



3 SPHERICAL ORBITS AND FRAME 
DRAGGING 

We now confine ourselves to the study of those orbits with 
constant r which are arbitrarily (not infinitesimally) lifted 
over the equatorial plane, i.e. with a finite value of Q. 

The conditions for the stability of a spherical orbit with 
radius r — ro are (see eq. (^|)) 

R(ro) = (9) 
f) Ft 
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Conditions (J^) and ( |lo| ) introduce two relations between 
r and the constants of motion E, L and Q, reducing the free 
parameters that characterize the orbit to two; thus, given a 
specific Kerr black hole (i.e. given the values of M and a), 
a spherical orbit is completely determinated, for example, 
by specifying its radius and the value of Q, whic h fixe s the 
amplitude of motion in the 9 direction (Wilkins (1972)). 



The motion is open, since the two fundamental frequen- 
cies v,f, and vg (proper or coordinate) are incommensurate; 
then the Fourier spectra of every function of the position 
of the test particle will contain a superposition of the two 
fundamental frquencies and of all their harmonics and will 
be of the kind 
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where (5 is an arbitrary phase. 

Therefore, the most natural signals to look for in such 
a system are the two fundamental coordinate frequencies 
themselves and the difference between them, that, as we 
will show, coincide with the unique correct definition of pre- 
cession frequency of the nodes of a spherical orbit. 

In fact we can compute exactly the coordinate period 
of the 9 motion: if we call 9± (with 9- < 9+) the two roots 
of the equation Q(9) = 0, we see from (^) that the particle 
oscillates on the coordinate sphere between the angles 9- 
and tt/2 + 9~. Dividing (|^) by (^) and integrating, we obtain 
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where /3 = 1-E 2 , k 2 



- / z+ (with z± — cos 2 9±) and K(k) 



and E(k) are the elliptic integrals of the first and second 
kind, respectively. 

The change of azimuth during one quarter oscillation of 
latitude is given by 



A<)> : 
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where n(fc) is the elliptic integral of the third kind. 

An orbit is called co-revolving (or prograde) if A0 > 0, 
counter-revolving (or retrograde) if Atfi < 0. If the 9 and cj> 
frequencies were the same, A(j> would equal it/2; it means 
that we can define 



Vg 2 



(14) 
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The angle by which the nodes of a spherical orbit are dragged 
during each nodal period is therefore 



AQ = 2tt\ 



(15) 



and, consequently, the coordinate precession frequency of 
the nodes (or Frame-Dragging frequency) is 



(16) 



An 

vfd = = \va — voi- 
re 

We stress here that this definition is different from the one 



given in eq. (2) of Cui, Zhang and Chen (1998) (vpd = 
v^AQ./2'k), with which it coincides only far from the source, 
where v^ ~ vq. But in this case it would be sufficient to 
consider t he weak-field limit, the w ell known Lense-Thirring 
equation (Lense & Thirring 1918), and this point clearly 
frustrates the aim of their work. In fact, approaching the 
horizon (to which the innermost stable orbit tends in the 
limit a — > 1) their definition leads to a divergence of the 
frame dragging frequency, while it can be shown that 

lim 2-kvfd = nsH 
where nBH is the angular velocity of the Black Hole 



(Christodoulou fc Ruffini 1971 Misner, Thorne fc Wheeler 



0. 



1973j) i.|e. the angular velocity of the Zero Angular Momun- 
turn Observers on the horizon. This is a relevant difference 
between this work and that of Cui et al. (1998). 

We chose to set, in our calculations, M = 7Mq, in or- 
der to compare our results with the 300-ffz QPO observed 
from GRO J 1655-40 (Remillard et al. 1997), the BHC for 
which the mass is most accurately measured (Orosz & Bai- 
lyn 1997). We considered only direct (i.e. prograde) orbits. 
In Figure 1 we plot the three frequencies calculated at se- 
lected radii as functions of a for Q — 1. The radii are the 
radius of the innermost stable circular orbit, which has been 
calculated solving the quartic equation 

d 2 R(r) 
dr 2 

r pea k — ri/rj (with rj slowly varying from 0.62 to 0.76 as a 
goes from —1 to 1), the radius of maximum surface emis- 
sivity of the disk (Page & Thorne 1974), 2r\ and 2r pea k- 
The value of Q = 1 (which corresponds to an 'opening 
angle' of the orbit over the equatorial plane which varies 
from about 3° for a = 0.5 to about 5° for a = 0.99) was 
chosen for simplicity, because, exploring the whole range 
Q — 0.01 — 10, we found the relative changes in the fre- 
quencies raising from ~ 2% (a = 0.5) to a maximum of only 
~ 5% (a — 0.99). It is immediately seen that for decreasing 
radii and increasing values of a the splitting of v$ and vg 
increases dramatically. Correspondingly the frame dragging 
frequency increases, reaching values that are comparable to 
ve for r = 2r pea k and 2n or even larger than ve for r = r pea k 
and ri. 

Figure 2 shows the three frequencies v^, ve and vfd 
versus v^, for selected values of the angular momentum of 
the black hole (a = 0.5, 0.9, 0.95, 0.998). These graphs 
represent the frequency changes that would take place if, 
for a given black hole, the orbital radius of the precessing 
matter changed (see Section 4). 

These results, obtained in a fully general-relativistic 
framework, admit a simple interpretation in terms of a New- 
tonian analogy. In the classical gravitational potential due to 



a spherical star, oc 1/r, the frequencies of motion for a bound 
orbit in the azimuthal (<j>), radial and latitudinal (9) direc- 
tions are all equal; this well-known property assures that all 
orbits close in this potential. Whenever a small perturbation 
is introduced, such as that due to the a star's oblateness, 
this property is lost and the ^-frequency ve becomes differ- 
ent from the (^-frequency v^. Then the spectrum emitted by 
a source on this orbit will contain all harmonics of the type 
nvtf, + mve, with n, m integers; of these, the line with, most 
likely, the largest amplitude is that at frequency ve — v^. 
This, in particular, is the classical precession frequency due 
to a Newtonian star not being perfectly spherical. Since in 
classical mechanics departures from spherical symmetry are 
always modest, we always find v$ — ve <C ve,v^. But in 
the gravitational field around a fastly rotating black hole 
such departures are much more significant, implying that 
this inequality is no longer satisfied. In other words, as de- 
partures from a Newtonian potential increase, whether be- 
cause we are moving to a strong-field limit, or because the 
black hole is rotating faster, we expect to move toward a 
situation where v$ — ve ve * v#. This is exactly what we 
see happening in Figures 1 and 2. 



4 APPLICATION TO BLACK HOLE 
CANDIDATES 

By using the black hole mass and angular momentum 
which have been measured (or indirectly inferred) for several 
BHCs, Cui et al. (1998) find a reasonably good agreement 
of the predicted point-mass precession frequencies at the ra- 
dius where the disk emissivity is highest, with the observed 
QPO frequencies. If disk precession is not confined to such a 
radius and differential precession takes place at a frequency 
close to the local frame dragging frequency, it remains to be 
demonstrated that a sufficiently narrow QPO peak matching 
the observations can be generated as a result of the differ- 
ent precession frequencies that might take place at different 
radii. More crucially, it is well-known that in viscous ac- 
cretion disks Lense-Thirring precession of the whole disk is 
strongly damped (Bardeen and Petterson 1975, Pringle 1992 
and references therein) . However there appear to be preces- 
sion modes that are strongly confined to the innermost disk 
regions and only weakly damped (Markovic & Lamb 1998); 
these are currently being investigated in greater detail. The 
mechanism reponsible for the excitation of these modes re- 
mains an open question. The perspectives for some kind of 
resonant excitation driven from an azimuthal asymmetry do 
not appear promising in consideration of the black hole "no 
hair theorem". 

An alternative possibility is that in the innermost disk 
region there are individual blobs moving like test-particles 
in the BHC field, executing also LT-precession, and mod- 
ulating the observed X-ray flux either through occultation 
or because they are self-luminous. The existence of discrete 
blobs is of course not an embarrassment for this argument, 
since their existence is required in all scenarios trying to ex- 
plain QPOs, in particular those involving weakly magnetic 
neutron stars in Low Mass X-ray Binaries (LMXRBs). This 
tendency of the disk to form discrete blobs seems to be inde- 
pendent of the nature of the accreting source; for instance, 
Krolik (1998) suggests that it be due to a need to circum- 
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Figure 1. The three coordinate frequencies (i/^, dashed line, i/g, dot-dashed line and vpn, solid line) of spherical orbits at four diffrcnt 
radii as functions of the dimensionless angular moment of an M = 7Mq black hole, r, is the radius of the innermost stable orbit and 
r P eak is the radius of maximum surface emissivity of the disk. Q is always set equal to 1. 



vent the local instabilities of the disk, irrespective of the 
properties of the accreting source. So this tendency might 
be present in accretion disks surrounding both neutron stars 
and black holes. 

By pushing further the analogy with neutron star 
LMXRBs, where high frequency QPOs are very often ob- 
served and succesfully interpreted in terms of the Keplerian 
frequency (note that in LMXRBs v$ ~ vg, see below), one 
would conclude that if blobs' precession is responsible for 
the QPO observed in BHCs, there is no obvious reason why 
QPOs reflecting the cf> and ^-components of the orbital mo- 
tion should not be there. Indeed, if, according to the model 
of Cui et al. (1998), the ~ 300 Hz QPOs of GRS J1655-40 
originate from the frame dragging frequency of blobs off the 
equatorial plane in the innermost disk regions, then the <j> 
and 9 frequencies of the orbital motion are v$ ~ 970 and 
950 Hz and vg ~ 670 and 650 Hz, while a ~ 0.88 and 
0.96, respectively in the case in which frame dragging QPOs 
are produced at the innermost stable orbit or the radius of 
highest disk emissivity. The difference between v$ and vg is 
large and two well-separated QPO peaks might be expected. 
These signals, however, have not been detected yet. Similar 



considerations would apply to the case of the ~ 67 Hz QPOs 
from GRS 1915+1050 (Morgan et al. 1997). 

The application of beat frequency models, BFMs, to 
those neutron stars systems that show twin kHz QPO 
peaks, allows us to identify the higher frequency kHz QPO 
(~ 800 — 1200 Hz) as arising directly from the Keplerian 
motion of blobs at the inner egde of the disk (moreover the 
neutron star spin frequency is inferred from the difference 
frequency of the twin kHz QPOs). Stella & Vietri (1998) 
noticed that the precession frequency of these blobs, as de- 
rived from the neutron star parameters inferred from BFMs, 
agrees well with a broad peak around 20 — 35 Hz that is ap- 
parent in the power spectra of three sources. In this model, 
therefore, v^—vg ~ 20—35 Hz, a separation that is compable 
to or smaller than the width of the higher frequency kHz 
QPO peak. Therefore it is not surprising that the signals at 
v$ and vg are difficult to disentangle in the case of neutron 



t The harmonic content of the three fundamental frequencies in 
the problem at hand will depend on the mechanism responsible 
for the generation of the signal(s) (e.g. self-luminous vs. occulting 
blobs) and on geometry, and is beyond the scope of this Letter. 
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Figure 2. The three coordinate frequencies {y$, dashed line, vg, dot-dashed line and vpr), solid line) of spherical orbits with Q = 1 for 
selected values of the angular momentum of the black hole (a = 0.5, 0.9, 0.95, 0.998) 



star LMXRBs. It should also be noticed that, around neu- 
tron stars with weak magnetic fields, a natural mechanism 
exists to lift matter off the equatorial plane, through the 
interaction with a spinning, tilted magnetic dipole moment 
(Vietri and Stella 1998). 

In the frame dragging interpretation of BHC QPOs, 
very specific predictions are also made in relation to the 
changes in v$ and ve that results from changes in the frame 
dragging frequency v$ — v$ (cf. Fig. 2). This would provide 
indeed a sensitive diagnostic to confirm the interpretation 
and study with unprecedented detail the motion of matter 
close to event horizon of a Kerr black hole. 

The work of L.S. is partially supported through ASI 
grants. 
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